An essential mechanism of heat dissipation in carbon nanotube electronics 
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Excess heat generated in integrated circuits is one of the major problems of modern electronics. 
Surface phonon-polariton scattering is shown here to be the dominant mechanism for hot charge 
carrier energy dissipation in a nanotube device fabricated on a polar substrate, such as SiC>2. Using 
microscopic quantum models the Joule losses were calculated for the various energy dissipation 
channels as a function of the electric field, doping, and temperature. The polariton mechanism 
must be taken into account to obtain an accurate estimate of the effective thermal coupling of the 
non-suspended nanotube to the substrate, which was found to be 0.1-0.2 W/m.K even in the absence 
of the bare phononic thermal coupling. 



Introduction 

Nowadays the vast majority of human activities are 
supported by advanced information technology which 
cannot become a transformative power without progress 
in semiconductor electronics. Exponentially increasing 
dissipated power density in integrated electronics circuits 
is one of the "grand challenges" of modern electronics 
[TJ |2] . Excess heat generated by electric currents in the 
elements of circuits cannot be completely removed be- 
cause of the high power density, which leads elevated de- 
vice operation temperatures, performance reduction, and 
ultimately results in hardware failures. 

Heat dissipation takes place via the transfer of thermal 
energy that is generated by hot charge carriers, to the 
cooling units eventually. The three classical mechanisms 
of heat transfer are convection, conduction, and radia- 
tion. The former is inapplicable in existing circuit archi- 
tecture having no fluid or gas flowing inside. The latter 
has a total emissive power proportional to T 4 according 
to the Stefan-Boltzmann law, and is not very effective 
within the acceptable range of temperatures [3 H] . Thus 
the major mechanism responsible for the device cooling 
is the thermal conductance. High thermal conductance 
can be achieved in bulk solids, but even a small vacuum 
gap is a quite effective thermal insulator. 

This is the case for transistor channels involving a 
nanostructure, for example a nanowire or a nanotube [5, 
H3|7], which is not chemically connected to the substrate, 
instead it is bound by relatively weaker van der Waals 
forces. This weak binding has a positive effect on the 
mobility of the charge carriers, but also has a negative 
effect on the thermal conductance. In existing nanotube 
transistors only the contacts have a good thermal ex- 
change rate while the thermal conductance over the van 
der Waals separation gap is very low, of the order of 
0.05-0.2 W/K.m EHl HD H2] • In addition, the to- 
tal thermal conductance of the contacts is not too high, 



being proportional to the cross-sectional area of the nan- 
otube. 

In this work we propose that a fourth mechanism of 
heat dissipation from the hot charge carriers in the nan- 
otube channel into a polar S1O2 substrate should operate 
and argue that such a mechanism would dominate the 
thermal transfer especially at large device currents (high 
bias voltages). Understanding the importance of such a 
mechanism invokes concepts that are specific to carbon 
nanotube (NT) field-effect transistors (FET): low ther- 
mal conductance due to the van der Waals separation 
gap (and possible surfactant coating) and the small con- 
tact area, high Fermi velocity of the NT charge carriers 
and the existence of the special electromagnetic (EM) 
surface modes in the vicinity of the NT channel. Hot 
electrons, as we will discuss below, can excite surface 
polariton modes. This mechanism while important for 
nanoscale channels is ineffective in a macroscopic system 
because the EM field of the polariton decays exponen- 
tially with the distance from the surface plane |T3 . This 
leads to a power low scaling of the scattering rate with 
the distance from the surface plane [T3J. A surface po- 
lariton mode would quickly dissipate into bulk substrate 
modes, thus providing a heat transfer mechanism with an 
efficiency exceeding the standard thermal conductance or 
radiation mechanisms. 

The existence of surface polariton modes was predicted 
in 1899 [TS] and were later observed at AM radio-wave 
frequencies. In condensed matter, a metal surface is 
known to support surface EM modes that is plasmon- 
polaritons; the focus of recent plasmonics studies |16j . 
The plasmon-polariton EM mode exists if the real part 
of the dielectric function e of the substrate equals nega- 
tive one |17j . Similarly, one can observe a surface polari- 
ton in a polar insulator where a strong optical phonon 
mode provides the required condition: e = — 1. Then 
the frequency of the surface phonon-polariton (SPP) is 
close to the frequency of the optical phonon. For exam- 
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pie, for a SiC>2 surface the optical phonon modes that 
can support SPP are at ZiwM = 50, 62, 100 and 149 
meV [T8j [El [20]. Hot charge carriers in a NT FET with 
energy exceeding the polariton energy would effectively 
release energy into the SPP channel. 

We stress that despite the fact that this mechanism 
is similar to a classical one for hot charge carrier en- 
ergy relaxation whereby the energy dissipates into opti- 
cal phonons of the NT EU El El EI], the SPP mech- 
anism transfers the thermal energy directly into a bulk 
substrate, while the classical one only redistributes it be- 
tween the NT electrons and the NT lattice, to be fol- 
lowed by thermal conduction into the leads and/or the 
substrate. The SPP mechanism also differs from the clas- 
sical radiation since the charge carriers in the channel 
couple to the evanescent EM modes only. The greatly 
enhanced surface electric field of such modes has much 
larger overlap with the NT channel than it would be in 
the case of the free vacuum EM modes [4]. 



Simulation approach 

We have adopted the model described in detail in 
the previous works [TU EH EH]- In brief, a standard 
tight-binding Hamiltonian for the NT charge carriers, 
a\a,j, is supplemented with Su-Schrieffer- 
Heeger terms for the electron-NT-lattice-phonon cou- 
pling (Eq.(l) of Ref-Ei): 



H NT = J2 K q 



+ a qi i \ [c v k+q c v ^ c c k+q c c ^ k 



here the matrix element of the interaction is M kq oc 3 = 
5.3 eV/A. The coupling constant is S oc t — t a , where 
t a = 3 eV is the bare tight-binding hopping integral be- 
ing modulated by the phonon modes. a qfl and c c / vk are 
annihilation operators for NT phonons (/it is a phonon 
branch index) and charge carriers in conduction and va- 
lence bands, respectively, that are labeled with the one- 
dimensional axial momentum q or k. In addition to the 
NT lattice phonons we include the operator of the inter- 
action with the SPP mode (cf. Eq.(2) in Ref.[27]l: 
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where a q ^ q±M is the annihilation operator for the z/-th 
surface phonon mode (not to be confused with the NT 
lattice phonons) , and an explicit form of V kq q± was de- 
rived in Ref . [27] : 
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here q and q± are the components of the SPP momentum 
that are parallel and normal to the NT axis respectively 



(see FigJTJ;) , m is the angular momentum transfer to the 
SPP mode (equals in case of the NT intra-subband 
scattering which is the strongest scattering channel), R 
is the NT radius, h ~ R + A A is the distance from the 
NT axis to the surface, I m (x) is the Bessel function of 
imaginary argument |35] . e is the elementary charge and 
T v characterizes the strength the electric field of the SPP 
mode [T7] : 
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where S is the normalization surface area, id^o ^ s ^ nc 
frequency of the longitudinal optical phonon of the polar 
substrate with the static permittivity eo and the low- 
frequency dielectric function e^. We should emphasize 
that while there are four surface modes present (the high- 
est energy one being doubly degenerate), only two of 
them have strong coupling strength [T3JE7]: SF 2 =1 5 = 
0.042,0.38, 0.069, 1.08, 1.08. The characteristic distribu- 
tion of the electric field of the SPP mode above the sub- 
strate surface is shown in Fig[l}i. The field oscillates 
along the surface (along q) and decays exponentially in 
the space (along The transverse component of the 
SPP momentum, qj_, is not conserved and has to be in- 
tegrated for the final result. 

To simulate the current-voltage curves we solved nu- 
merically the steady-state Boltzmann equation for the 
NT in a constant electric field (including as many sub- 
bands as needed for given temperature and drain field 
values). One can assume that in the bulk of the FET 
channel the drain voltage Vd is related to the applied elec- 
tric field via the effective channel length L: F — Vd/L. 
The scattering rates for the Boltzmann equation are as 
follows (see Refs.[14] and [27] for further details): 
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here fiuj^Q is the frequency of the SPP mode [T7J; n qv 
its phonon occupation number, Efu(k) is the energy of 
the final/initial state as obtained from the Hamiltonian 
H e . Non-equilibrium distribution function g k has been 
calculated numerically for every given value of the tem- 
perature, field and charge density and then it was used 
to compute the total electron current, Id, as well as to 
determine the partial phonon emission rates for the heat 
dissipated in the NT lattice Pj and that dissipated di- 
rectly into the substrate via the SPP mode P — Pj, ac- 
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FIG. 1: Role of SPP inelastic scattering mechanism for 
NT FET transport. (a) Current-electric field curves for 
[17,0] NT at the doping level p = 0.1 e/nm, T = 
77,150,210,300,370,450 K as indicated by the color code. 
Full line corresponds to the room temperature characteristic 
with both SPP and NT phonon scattering mechanisms in- 
cluded, as compared to the SPP scattering only (red dotted 
curve slightly above) and NT phonon scattering only (light 
blue dotted curve at the top). Insets in (a) shows the drift 
velocity vs. the doping level at F — 5 V/p,m and T = 300 K. 
(b) The energy relaxation time vs. applied electric field us- 
ing both SPP and NT scattering for the same NT and same 
temperatures. Inset in (b) shows dependence of the energy 
relaxation time on the doping level at F = 5 V/p,m and 
T = 300 K. (c) Schematics of the NT channel on a Si0 2 
substrate (3D image is generated using Molecular Dynamics 
in NAMD/VMD). Accepted convention for the components of 
the in-plane wave- vector of the SPP mode and the coordinate 
system are shown, (d) Distribution of the electric potential 
and its gradient for a Surface Phonon Polariton (SPP) mode. 
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where P = I^F are the total Joule losses. A similar 
equation holds for NT losses Pj with the scattering rate 
obtained from Eq. (JTJ. In the following we study depen- 
dencies of the SPP to the NT channel loss ratio, defined 
as £ = P/Pj — 1, on the bias, density, and temperature. 



SPP Scattering and NT-FET High Bias Regime 

The role of the SPP modes with respect to charge car- 
rier scattering (so-called Remote Interface Phonon scat- 
tering) for Si devices has been first studied by Hess and 
Vogl |25j , though this mechanism was not found to dom- 
inate the transport [26] . As we show in this work, for 
nanotube devices the SPP scattering mechanism is very 
fast, as been indicated already in our earlier paper |27j 
and cannot be neglected even for the low-field transport 

El- 

The SPP energy relaxation rate exceeds the intrinsic 
NT phonon relaxation rates (for both optical and acous- 
tic modes) for the whole range of the charge densities 
and applied electric fields (gate and drain biases) studied 
here as shown by the results in Fig. [T]d. 

Recently SPP scattering was demonstrated to deter- 
mine the mobility in graphene on a polar substrate and 
to limit the ultimate performance of graphene devices 
[2^1 I3U1 |3~T] . FigjTJi shows typical current-electric field 
curves calculated with and without the SPP mechanism 
taken into account. Neglecting the SPP channel sig- 
nificantly overestimates the low-field current. Two im- 
portant conclusions are drawn from these data: firstly, 
the negative differential resistance region of the current- 
voltage characteristic is totally removed when the SPP 
channel included, because the very effective SPP energy 
relaxation prevents the hot electron run-away as we fur- 
ther discuss in Ref. [TJ. Secondly, the SPP mechanism 
dominates over the NT phonon scattering for both low- 
and high-field regimes but the remnants of the non-SPP 
behavior can be seen at low applied field especially at low 
temperature. 

Results on the relative importance of the SPP and NT- 
phonon scattering channels for the energy relaxation of 
the hot charge carriers are shown in Figj2] In Fig(2^a) 
we plot, in a logarithmic scale, the partial specific en- 
ergy losses due to the SPP channel (full curves) and only 
the NT phonon channel (broken curves), calculated for 
[17,0] NT at the doping level p = 0.1 e/nm at various 
temperatures. The SPP losses are up to two orders of 
magnitude larger than the NT phonon losses. Figj2^b) 
presents the ratio of the SPP partial losses to the NT 
phonon partial losses calculated for the same NT. The 
high frequency NT optical phonons increase the NT par- 
tial losses at higher fields and higher T (Figpb). Even 
in that limit their ratio £ is always greater than one. 

We also studied the dependence of the loss ratio on the 
doping level p (Figj2^, inset) at the fixed applied electric 
field F — 5 V/p,m. The relative decrease of the SPP con- 
tribution at higher electron density (higher gate bias) is 
shown in Figj2^e) as a percentage of the total losses and 
is due to the SPP and NT phonon losses having a dif- 
ferent functional dependence on p. The non-equilibrium 
distribution function (obtained as a numerical solution 
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of the Boltzmann equation) varies with the doping level. 
The SPP losses, correlated with the non-equilibrium dis- 
tribution function, scale similar to the total losses (pro- 
portional to the current) and are approximately a lin- 
ear function of the doping level (total number of elec- 
trons). The NT phonon losses are superlincar in p and 
grow faster, thus their ratio £ decreases with p. 

To understand the dependencies of the NT and SPP 
losses on field strength and temperature, we show the 
scattering rate of electrons via these two channels in dif- 
ferent bands as a function of energy (see Supplemental 
Information Fig. 3). In the first two energy bands the 
SPP scattering dominates over the NT scattering, while 
in the higher energy bands the NT scattering becomes 
comparable to that of the SPP scattering. This is because 
in the SPP scattering the angular momentum is con- 
served Am = 0, so that only intra-band scattering takes 



place, while NT optical phonons can lead to both intra 
and inter-band scattering and the phase space for such 
scattering grows with the band index. Therefore, as the 
electronic distribution becomes hotter either due to the 
field or temperature the relative role of the NT phonons 
increases and the loss ratio £ decreases. Fig. 4S (Supple- 
mental Information) further supports that the electronic 
distribution function rather than the scattering rate (see 
Eq. ([6|) is primarily responsible for the value of the loss 
ratio at the high bias regime. Therefore the NT and SPP 
losses scale in exactly the opposite way with temperature. 



Near-Field Thermal Conductance and Self-Heating 
of NT-FET Channel 



Next we discuss another aspect of SPP scattering, 
namely we address how the heat generation is controlled 
by the fast SPP energy relaxation rates. 

So far we set the NT temperature equal the substrate 
temperature. That is we assumed an ideal thermal cou- 
pling to the substrate which is a condition not confirmed 
experimentally. As a next level of approximation we 
use thermal coupling values (per NT length) from the 
literature HOj HU [12] ranging from g a = 0.05 to 
0.2 W/K.m and estimate the FET channel heating. The 
lower bound of the thermal coupling, as it will be shown 
below, gives so small a thermal exchange to the substrate 
that it is almost equivalent to assuming no thermal cou- 
pling at all (see Supplementary Figure 1). Thus we an- 
alyze first the temperature of the NT FET channel for 
the high value of the thermal coupling to the substrate 
g Q = 0.19 W/{K.m). 

All energy loss of the hot electrons, same as the total 
dissipated power, is distributed between two channels: 
the SPP losses which are to be subtracted from the total 
heat flux in the NT, and the NT losses which have to be 
partially transmitted to the substrate via the coupling 
g . Fig|2jc) presents an effective channel temperature 
(the Joule overheating) as a function of the applied elec- 
tric field calculated for NT phonons only (upper curves' 
family) and for both NT and SPP channels (lower curves' 
family). The SPP scattering channel couples to the sub- 
strate phonons and thus transfers most of the excess ther- 
mal energy of the hot electron directly into the substrate. 
In contrast, the NT phonon scattering transfers the en- 
ergy to the NT lattice. As a result the steady-state tem- 
perature of the NT lattice increases from the ambient 
temperature, T S11 &. For a very long NT FET channel all 
edge effects due to the electron and phonon thermal con- 
ductivity can be neglected, as well as a heat flux due to 



the hot electron current into the drain electrode. Given 
the thermal coupling rate we estimate the steady-state 
temperature as: 

T = T sub + AT = T sub + -± (7) 

9o 

where Pj is the dissipated power per NT length due to 
the NT phonon scattering. The Joule overheating AT 
is found to be orders of magnitude smaller compared to 
the case when SPP channel is neglected and all the mea- 
sured Joule losses are attributed to the NT phonons, i. 
e. Pj = P = I d F (see Fig(2j:). The ratio of the SPP 
losses to NT losses, £, as a function of the applied field 
and doping level is given in Figj2j) and in the inset (e). 
The temperature rise AT is inversely proportional to g Q 
according to Eq.([7]). Even though the loss ratio depends 
on g Q (see below), which is not precisely known, it is 
safe to conclude that the temperature of non-suspended 
NT channel will be much lower for the same total power 
density P. 

Within this non-self-consistent scheme one overesti- 
mates the overheating because one neglects AT when 
calculating the current, the SPP and Joule heat (NT) 
losses. Next we part this approximation and provide fully 
self-consistent analysis of the channel temperature. 

Given the expected temperature rise of the NT FET 
channel (for a given value of the bare thermal coupling) 
we recalculate the electron and phonon distribution and 
obtain the self-consistent NT temperature by iteration. 
At rt-th step we solve the Boltzmann equation using the 
NT-lattice temperature T„_i, calculated at the previous 
step according to Eq. ^ , and then compute the total en- 
ergy released in the NT-phonon subsystem. Next, using 
Eq.(|7]) we calculate the new lattice temperature T n . The 
process is iterated until convergence. 
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FIG. 2: Comparison of the hot electron energy relaxation via 
SPP and other inelastic scattering channels, (a) Logarithm of 
the partial specific energy losses due to the SPP channel (full 
lines) and the NT phonon channel (broken curves) for [17,0] 
NT at the doping level p — 0.1 e/nm (the temperature range 
is indicated by the color code T = 77, 150, 210, 300, 370, 450 K 
arrows point to increasing temperatures), (b) The ratio of 
the SPP partial losses to the NT partial losses, £ = P/Pj — 
1, vs. applied electric field for the same NT. (c) Estimated 
effective NT FET channel temperature vs. applied electric 
field calculated for [17,0] NT on Si0 2 substrate. NT phonon 
contribution (upper curves' family) is compared with both NT 
and SPP channels (lower curves' family) for the bare thermal 
coupling to the substrate g = 0.19 -g^r, p = 0.1 e/nm. 
(d) The loss ratio £, calculated including self-consistent NT 
heating, is shown as a function of the doping level and the 
applied electric field. Here g — 0.19 g^r, T BU t = 300 K. (e) 
Relative contributions of the SPP (red curve) and NT (blue 
curve) channels to the total energy losses (black curve) vs. 
the doping level at F = 5 V/ fim; the inset shows their ratio. 



At zero bare thermal coupling go — 0, SPP chan- 
nel controls the Joule heating. The field and concen- 
tration dependence of the FET temperature at g = 
is presented in Fig|3^. The upper (purple) curve corre- 
sponds to [17,0] NT on Si0 2 substrate at T sub = 300 K, 
p = 0.4 e/nm. A temperature rise higher than 600 K 
can be expected at such conditions [32 . By decreas- 
ing the doping level 4 times (red curve) we decrease the 
steady-state temperature due to a smaller value of the 
Joule losses, proportional to the current ~ Id ~ p. Addi- 
tional decrease is associated with a gradual dependence 
of the loss ratio £ on the density (Fig{2^ inset). To fur- 
ther explore the nature of this cooling down of the NT 
FET channel we compute the effective thermal coupling 
to the Si0 2 substrate, defined as g = (I d F)/AT. We 
extract it from the self-consistent temperatures data (as 
in Fig j3^i) and plot it vs. the applied field, F, and vs. 
the total dissipated power per NT length, P, in Fig(S]3. 



In a real device the non-zero bare thermal coupling 
g a > results in even higher values of the total ther- 
mal coupling. Effective thermal coupling g and the ratio 
g/g Q are presented in Fig|3] panels (c-f) vs. the dop- 
ing level at fixed electric: field F — 5 V/p,m and also vs. 
the field for p — 0.1 e/nm and various bare couplings 
g a = 0.05,0.10,0.15 and 0.19^ (from red to blue). We 
conclude that the SPP, near-field, thermal conductance is 
always larger than the bare thermal coupling by at least 
an order of magnitude. 

We note that in the presence of dissipation mechanisms 
to both SPP and NT channels, the total thermal con- 
ductivity of the NT FET becomes a function of bias and 
density and the analysis of the electric heating through 
Joule's, Fourier's and Ohm's laws, often used in litera- 
ture [33l [34] , must account for these dependencies. In 
such a case, one has to calculate a self-consistent non- 
equilibrium distribution function first and find the lattice 
(phonon) temperature. The total thermal flux has two 
components: the SPP channel flux and the NT-substrate 
flux. We predict that for supported (non-suspended) NT 
device on a polar substrate the SPP dissipation chan- 
nel is major and dominates over the NT-substrate ther- 
mal dissipation. Therefore, an effective thermal coupling 
through the NT-substrate interface appears to be signif- 
icantly larger as compared to the bare thermal coupling 
on a non-polar substrate. 

Conclusions 

In this paper we presented a microscopic quantum 
modeling of a novel heat dissipation mechanism for nan- 
otube electronic devices. This mechanism is specific for 
NT devices fabricated on polar substrates, such as SiC-2, 
due to (i) existence of surface EM modes at the frequen- 
cies of the surface phonon-polaritons, (ii) strong coupling 
of such modes to the charge carriers in the NT lying 
on the substrate. We note that a similar SPP thermal 
coupling should exist in other ID and 2D systems fab- 
ricated with other channel materials on polar insulator 
substrates. In particular, high-k oxides are expected to 
produce strong SPP scattering due to their large Frohlich 
constants and their low optical phonon frequencies. 

Using a semiempirical quantum approach we have cal- 
culated current-voltage curves with and without SPP 
scattering and concluded that the SPP mechanism dom- 
inates the scattering and determines the drain current 
in the whole range of drain voltages, and for all stud- 
ied doping levels and temperatures. The current (charge 
density) and temperature scaling of the SPP and NT 
scattering mechanisms are different, thus allowing verifi- 
cation of our predictions experimentally. In this work we 
focused on the FET high bias regime and demonstrated 
that the SPP mechanism, being much more efficient than 
the scattering by NT lattice modes, results in the elim- 
ination of the negative differential resistance part of the 
current- voltage curve. We analyzed the relative impor- 
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FIG. 3: SPP as a heat dissipation mechanism, (a) Self-consistent effective NT FET channel temperature vs. applied electric 
field calculated for [17,0] NT on SiC>2 substrate, T 3u i = 300 K, the doping level is p = 0.4 e/nm (purple, upper) and 0.1 e/nm 
(red, lower). A bare thermal coupling is set to g = 0. (b) Logarithm of the effective thermal conductivity via the NT-SiC>2 
interface for the same NT as in (a) due to the SPP coupling, as extracted from the self-consistent temperatures data, vs. 
applied electric field. Inset shows the same data vs. dissipated power per NT length, (c-f) Effective thermal coupling g and 
g/g -ratio vs. the doping level at F = 5 V/p,m, (c) and (e); the same vs. applied electric field at p — 0.1 e/nm, (d) and (f). 
Bare thermal coupling is g a = 0.01, 0.05, 0.19 in (c) and (d) and g a = 0.01, 0.05, 0.10, 0.15 and 0.19 (from blue to red) 
in (d) and (f). 



tance of the SPP channel, comparing to all the other NT 
phonon mode scattering channels. In the whole studied 
range the SPP channel dominates. Thus, most of the en- 
ergy losses are dissipated directly into the polar substrate 
and do not contribute to the FET temperature rise (Joule 
overheating) . We showed that the SPP thermal coupling 
increases the effective thermal conductance over the in- 
terface between the NT and such polar insulator as S1O2 
by an order of magnitude. The dependence of the effec- 
tive thermal conductance of the NT/SiC>2 interface on 
the channel doping level was computed and may be used 
to verify the model predictions. We note that plasmon- 
polaritons, existing at a metallic surface, may also result 
in a strong thermal coupling to the nanotubes. 
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